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> ■ Abstract 



0{h) effects that modify the classical orbit of a charged particle are described 



■^ . for the case of a classical spin- 2 particle moving in a constant magnetic field, 

ON 

c/3 . using a manifestly covariant formalism reported previously. It is found that 

"^ . the coupling between the momentum and spin gives rise to a shift in the 

O , cyclotron frequency, which is explicitly calculated. In addition the orbit is 

a ■ 

found to exhibit 0{h) oscillations along the axis of the uniform static magnetic 



^ . field whenever the gyromagnetic ratio (7 of the particle is not 2. This oscillation 

is found to occur at a frequency oc | — 1, and is an observable source of electric 
dipole radiation. 
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I. INTRODUCTION 

Closed orbit correction is a critical component of background reduction in high energy 
electron and proton accelerator experimental detectors. In a collider ring, tight quality 
control of the closed orbit is always essential to the efficient operation of the detectors. 

Higher-order non-linear relativistic processes may measurably effect the dynamical tra- 
jectory of a classical spin-i electron in an electron accelerator or a proton in a particle beam. 
Such effects can in principle be seen after many orbits in a long-lived charged particle beam 
in a storage ring when conditions for resonance are met. For example, a key issue in collider 
physics is to what extent the interaction of non- vanishing magnetic field gradients along the 
particle's trajectory with the intrinsic magnetic dipole moment of the particle drives the 
particle from the ideal design orbit (the relativistic Stern-Gerlach effect). This may alter 
the depolarizing resonance strengths and widths, and may induce new resonant modes that 
effect the spin. Hence a quantitative description of the contribution of non-linear relativistic 
effects to closed orbit control is desirable. After over fifty years of study, this subject con- 
tinues to be an active research topic 0] - 0. Failure to study these processes could lead to 
a serious deficiency in our knowledge of the details of high-energy classical dynamics. For 
example, small innocuous non-linear terms can sometimes lead to spectacular observable 
experimental consequences (recall the hysteresis in cyclotron resonance based on a weak 
non- linear relativistic mass effect). 

Moreover the recent questions concerning the origin of proton spin (Ashman, et a/ [0, 
Close and Roberts [^, Ellis and Karliner P, Meng Ta-chung, et al |]lOl) and the observation 
of a strong spin-dependence of high-energy proton-proton interactions underline the need 
for a complete understanding of classical spin. 

In this paper we consider the classical motion of a charged particle with intrinsic spin 
orbiting in a circle about a constant magnetic field normal (in the zero^^ approximation) to 
the plane of the orbit. Such fields occur in an ideal accelerator or are approximated by the 
stellar magnetic field near a pulsar. We study the dynamical problem classically to 0{h) 



and derive several new experimentally testable predictions of classical dynamics. We also 
clarify a classic result, which states that the ratio of the spin precession angular frequency 
to the cyclotron frequency is given by 1 + 7(| — 1) Hi' o' B, where g is the gyromagnetic 
ratio of the particle and 7"^ = 1 — %. It is shown that this is a time- averaged result, 



and derive the time-dependent expression in Equation ||3.17|| . We shall adopt a classical 
description throughout this paper so that we may describe the electron with the same 
formalism as the proton. A quantum mechanical self-consistent description of an electron 
moving in the prescribed fields would require including radiative corrections to the electron- 
photon vertex in order to account for the fact that g y^ 2 for the electron. (Not so for 
a quantum mechanical description of protons, where one does not hesitate to introduce a 
phenomenological gyromagnetic ratio.) In the classical description, ^f is a parameter that 
one sets equal to an experimentally determined value. Although a classical description is 
many times inadequate, in this case we shall see that the classical description predicts and 
describes some effects that are difficult to calculate using relativistic quantum mechanics. 

Although not treated in this paper, the method employed to formulate this dynamical 
problem can be extended to handle the case of charged-particle beams. This is important 
since the actual synchrotron beam is comprised of multiple bunches of charged particles, 
with up to 10^^ particles per beam bunch. The dynamics of the beam is determined by 
the applied external fields plus the electromagnetic fields generated by the particles in the 
interacting bunches. The normalized precessional frequency is shifted by the beam-beam 
interaction and these shifts can also be approximated. 

II. EQUATIONS OF MOTION 

The Lorentz force equations for the four-velocity ^ = x° coupled with the Thomas- 
Frenkel-Bargmann-Michael-Telegdi (BMT) equations for the dynamical evolution of the 
Pauli-Lubanski spin vector govern the motion of a polarized particle in the limit h —>■ 0. 
Here r is the proper time. These equations have been extended to higher order in h. Al- 



though there are many specific models O" liJ, we shall employ the dynamical equations 
derived by Nash H . This coupled set of equations reduces exactly to the standard BMT 
and Lorentz force equations in the limit h —>■ 0. Moreover these equations are derived from a 
Lagrangian using a variational principle, which ensures that the dynamical equations do not 



conflict with conservation laws. We turn briefly to a summary of the approach of Ref ||19|| . 
Units are used in which the speed of light is one. The Minkowski spacetime metric is 
Va/3 = diag{l, 1, 1, —1). The dynamical variables are the particle's position x°', a = 1, . . . , 4 
in spacetime and ip, which is a real eight-component spinor. ip may be regarded as a column 
vector comprised of the direct sum of a real four- component Dirac spinor A [transforming 



under the real irreducible representation of 5*0(3, 3) defined by Dirac (see Ref [^)] and the 
transpose of another real four- component Dirac spinor ^ that transforms under the inverse 
irreducible representation 



V' 



'x\ 



(2.1) 



where the tilde denotes the transpose. ip defines an orthonormal tetrad in Minkowski space- 
time I21I. The members of the tetrad are constructed from a sum of products of the com- 
ponents of ■?/'. It is known c3 that ip may be regarded as a (split) octonian. 
Explicitly, the timelike member of the tetrad is given by 

E(l)(r) = -^^r^r>, (2.2) 

where the r„ matrices are real 8x8 analogs of Dirac's gamma matrices, and is parallel to 
the four-velocity x" when Equation ||2.4|| (below) is satisfied. In passing, we record for later 



use the definition of the third member of the tetrad 

i?f3)(r) = i^r^Fr"^. (2.3) 

For simplicity of presentation it is assumed that ip has been normalized to ipip = x^. The 
general case is obtained by dividing by a normalizing factor proportional to ^7^A. 



The classical set of dynamical variables {x°',ip} are not all independent. This is because 
the timelike E?^^ constructed from an arbitrary ip is in general not parallel to x". But in 
the free-field case one knows that only one timelike vector is required for the description of 
particle 'dynamics.' In order for {x°',tp} to describe a massive particle with spin, ip must 
be constrained so that E?^^ and x" are parallel. The crucial result we shall use is that the 
timelike member of the tetrad can be permanently aligned with the four-velocity x° of the 
particle by imposing the constraint H 

(Vax" + V-i;„x" r^) ^ = 0. (2.4) 

This classical constraint bears a striking resemblance to the quantum mechanical Dirac equa- 
tion. This suggests that one may quantize this theory by "quantizing the constraint" and 
postulating minimal coupling. However we do not treat the problem considered in Section 
PI quantum mechanically. To do so one must consistently take into account higher order 
corrections to the magnetic moment of the particle that arise from radiative corrections to 
the lepton-photon vertex operator. In the classical approach one may use the experimen- 
tally measured gyromagnetic ratio in the classical equations of motion. One obtains an 
approximation for the particle's orbit. 

The constraint Equation ||2.4|] is incorporated into the theory using a Lagrange multiplier. 
If this constraint is satisfied then the third member of the tetrad may be identified with the 
Pauli-Lubanski spin vector, and moreover, the intrinsic electric dipole moment of the particle 
vanishes in a rest frame T^apx^ = 0, where S"^ is the spin tensor of the particle defined by 

S"^ = --^f](M"^)^. (2.5) 

Here Q is the symplectic form on the spinor manifold and 

M°'^ = -i[r",r^]. (2.6) 

The Lagrangian for the theory is 



L = — Ma/— Xq-x" ipQip + eAax" + Lagrange multiplier term (2.7) 



where e is the charge of the particle, Aa is the electromagnetic vector potential, and M 
denotes the effective mass of the particle, which takes into account the spin-field interaction 
energy contribution to the particle's mass. It is given by 



M = m^l- £^S"^i^«. (2.8) 

where m is the rest mass of the particle (or simply a parameter with dimension of mass, 
which is to be renormalized to the observed rest mass) and Fap is the electromagnetic field 
tensor. Fajj contains a contribution from the radiation reaction of the particle. If this 
contribution to Fai3 is completely neglected then the equations of motion admit unphysical 
solutions in the free-field limit. 

The canonical four-momentum p^ is given by 

Pa, = Mxa + eA„ + hT.o.ph'^ (2.9) 

. the classical version The equations of motion derived from this Lagrangian are 

n^ = -hAF^"^ (i^^"/^ + ^"^/^) (2.10) 

where 

6" = x" - ^F%x^. (2.11) 

2M ^ ^ ' 

Also 

Mx° = eF^'pX^ + P"^T^, (2. 12) 

where P"'' is a projection operator P"'' = 7]'^^' + x'^x^' and T^ = f||S^^^ - ^ (fiL^pW^ 
contains the contribution to the force due to the Stern-Gerlach effect. 



For completeness we mention that upon substituting the result Equation ||2.1CI|| into the 
definition of the tetrad E?--. yields 

where j = 1, 2, 3 labels the spacelike members of the tetrad. It is worth noting that in the 
limit h -^ Equation p.l2| is the Lorentz force equation and Equation | |2.13| |, with j = 3, 



is exactly the BMT equation. 
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III. UNIFORM STATIC VERTICAL MAGNETIC FIELD 



In order to simplify the notation we define two sets of three 8x8 matrices 



a = {M^^M^\M'^) 



(3.1) 



and 



14 /I #24 .#34n 



m= {M^^,M'^,M 



(3.2) 



We employ the representation of the F matrices defined in References [Q and [^ . 

The electromagnetic field tensor Fajj is represented in terms of the electric field E and 
magnetic field B as 



aP 



(3.3) 



i?3 —B2 El 

-B, Bi E2 

B2 —Bi E2, 

—El — E2 — -c/3 

We assume that the h ^ Q limit of the Lorentz force equations have been solved so that 



v 



/ 



x" = x'^ij) is known. Substituting this result into the 



^ = V^(r)^, 



dh 



h^o 



limit of Equation |p.iq| yields 



(3.4) 



where 



^ = 2^[~2^"^ + ^f ~ 1)*'(*-^/^'^ - ^PP'^^)W^ 



is real 8x8 matrix that is a known function of tau. 



In the notation introduced above Equation |p.5|| reduces to 



-V 
e 



-{l+7'(|-l)k-B + (|-l){x-B}a-x 
- 7(- - l)m ■ B X X - (- - l){x ■ E} m ■ X 



.2 9 
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y (1 _ /?^^)m ■ E + 7(^ - l)a ■ X X E. 



(3.5) 



(3.6) 



For the limiting case h ^ 0, B = zB^, -B3 = constant, the electromagnetic field has only 
one non- vanishing independent component in the lab frame. This is F12 = B^. The Lorentz 
force equations ||2.12|| are independent of ip to zero^^ order in h. The zero^^ order solution 



/i ^ is of course well known. For a circular orbit in the xAy plane a solution to the Lorentz 
force equation is x^ = P'j cos^Uqt) , x^ = —P'ysm^UQr), x^ = and x^ = ■y = constant, 
where Uq = ^ , P = ^ = v = constant and 7^ = Y^ga- The cyclotron frequency is 

^cyclotron =\ ^0 \ ll- 

If we put T = ijJqt and v = 7(| — 1), and use Equation ||3.6|] , then Equation JS^ yields 



^ = -{(1 + yu)M'^ + /37z/[sin(T)Mi^ + cos(T)M24]}^. (3.7) 

Upon making the substitution 

= exp[(l + -fiy)TM^^]ip, (3.8) 

and simplifying, one finds that satisfies 

-^ = -/37i/[-sin(7i/T)M^^ + cos{-fiyT)M^'^]^ 

= -/37i.e^-^^^''M24e-^'^™'V, (3.9) 



where we have used 



= - sin(7z/T)M^^ + cos{-f iyT)M'^^. (3.10) 



The general solution to Equation |p.9|| is now easily found. One finds that ^jJ 



^B3,g,/3,r "0(0), where the constants of integration are in ip{0) and 
^^B3,5,/3,r — e e ' 

_ -K[cos(T)Mi4-sm(T)Af24] -(l+i/)TMi2 kM" /o i i \ 

and K = llnji^^}. e'^ is Bondi's 'k' factor Efl. ABs,g,(i,T is decomposed into the product of 
a boost in the direction of — x(0), times a r-dependent rotation about the guiding magnetic 



field, times another boost in the xcos(T) — y sin(T) oc x(T) direction. Ab3,c,,/3,t commutes 
with r^ and hence does not affect the polarization E3. 

In order to calculate the Pauli-Lubanski spin vector from ^(r) = 5*0^(0), with So = 
^B3,g,i3,T we first solve the eigenvalue problem Equation ||2.4|| for ip{0), choosing initial phases 
that simplify our work: 



r,x"(o) + v'-x«(o)x"(o) r^) ij{o) = o. 



(3.12) 

We find that there are four linearly independent solutions to this eigenvalue problem for 
x^ > 0, two solutions for spin up, and two solutions with spin down. As a check on 
our work we substitute ipi^) and Sq into the definition EZ){t) = —^ip{0)SoT^T°'Soip{0). 
Each distinct spinor solution, when substituted into this equation yields the four-velocity 
-£'(4) = {Plcos(T), — /?7sm(T), 0, 7) = x", as required. Next, substitution of the four spinor 



eigenvectors in turn into the definition Equation ||2.3|| yields two distinct Pauli-Lubanski 
four- vectors, which differ only in the sign of the third component. For example, we find that 
a "spin up" spinor eigenvector is 

i>{0)^ = (0, 0, 1, 0, 0, -/?7, 0, 7)/v^, (3.13) 



and ip = AB3,g,/3,T i'{0)^ = is given by 

/ 



^ 



-/57 sin(T/2) sin(z/T/2) 

-/57cos(T/2)sin(z/T/2) 

cos(T/2) cos{uT/2) - 7sin(T/2) sin(z/T/2) 

- sin(r/2) cos(z/T/2) - 7 cos(T/2) sin(z/T/2) 

-/57 sin(T/2) cos(z/T/2) 

-/37 cos(T/2) cos(z/T/2) 

cos(T/2) sin(z/T/2) + 7sin(T/2) cos{uT/2) 

- sin(T/2) sin(z/T/2) + 7 cos(T/2) cos(z/T/2) 



\ 



(3.14) 



The Pauli-Lubanski spin four-vector in this case is 



^fs) = (-2/^[(7 + 1) cos(z/ + l)r + (7 - 1) cos(z/ - l)r], 



-/3[(7 + 1) sin(z/ + 1)T - (7 - 1) sin(z/ - 1)T], 
7""*^, — /9^7cosz/T). 



(3.15) 



-E(^-) is just the Lorentz transform to the lab frame of the spin 3- vector s = (— /3cos[(z/ + 



l)T],/9sin[(z/ + 1)T], -), which may be verified by transforming E?^^ to an instantaneous 



rest frame with the Lorentz boost Lhoost defined as usual by x'* = x'^{6j — rfrij + rfrij) = 
{6'j — n^nj)x^ + n*7(n • x — Pet) = Lhoost *Q,a;° and ct' = •y^ct — /5n ■ x) = Lhoost a^;", where n' 
is a unit vector parallel to the 3-velocity of the particle. Applying this boost to E?'^^ yields 



j-irest frame a 
^ (3) 



-/3cos[(z/ + l)T],/5sin[(z/+l]T),-,0). 
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(3.16) 



Since XaE?^^ = 0, the angle Or between the three-velocity x^ and the three-spin Ej^^ is uT. 

def 

One sees that Or = vT = uujqt = 7(^2 



'£ — l)£^A which is the well known result for the 
precession of the longitudinal polarization c3l . 

The angular velocity of precession ujprecess that is measured in the lab frame is given by 



^ precess I ^cyclotron 



_, dE,n\ 



(=^) dT (^) dT 

^/2-l 
/52sin2(i/T)' 



E%.'^-^)/{El,,f + {E^^- 



\z)) 



(3.17) 



The average over a time vT = 27r is 



< ^precess / ^ cyclotron >— 1 + ly— — 1) — 1 + I^. 



(3.18) 



which is the well-known expression often identified as the normalized spin precessional fre- 
quency (see, for example, References ^, and 0). 

In order to write out and solve the momentum equations [ 2.12|| to 0{h) we must first 
evaluate the components of the spin tensor. Substituting ip form Equation [ ^.14| ] into 
Equation [ Ol ] yields 
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^.3\ / 



5.14 
5.24 
5.34 



V^ / 



P[— cos(T) cos(z/T) + 7 sin(T) sin(z/T)] 

/3[sin(T) cos(z/T) + 7cos(T) sin(z/T)] 

1 

-/3sin(T) 

-/5cos(T) 

-/527sin(z/T) 



\ 



(3.19) 



In terms of their Fourier decompositions, S^^ = |[— (7+I) cos((z/+l)T) — (7 — 1) cos((z/— 1)T) 
and S^i = f [(7 + 1) sin((z/ + 1)T) + (7 - 1) sin((z/ - 1)T). 

A word about normalization. Throughout this paper the normahzation of ip is deter- 
mined by the arbitrary requirement that the i^f^x comprise a normalized set of four mutually 
orthogonal vectors, which is reflected in the result that S^^ = 1. This is not a statement 
about the magnitude of the classical spin. The magnitude of the classical spin is not pre- 
dicted by this theory. Instead it must be imposed as an initial condition; in virtue of the 
equations of motion this magnitude is a constant of the motion. Given this magnitude and 
g we may employ this formalism to compute the trajectory. Conversely, one can impose 
the observed spin magnitude as an initial condition and then use the predictions of this 
theory and experimental measurements to determine g. One sees that to correctly apply 
this dynamical formalism to a spin-s particle one must arrange that S^^ = s. This is easy 
to accomplish by imposing this as an initial condition, and is manifested in the simple 
replacement ip{0) 1-^ ^ip{0), which we shall henceforth employ. 

We turn now to the 0{h) solution of the extended Lorentz equations for the momentum. 

m — s^hujQ to 0{h). Also to this order, 



The effective mass M 



1 - I^S"'^ P 



mx/i-^^z.--rat3 



hb"" = —GUujqF^'p qx^ /B3, where G = | — 1 and qX^^ refers to the 0(1 
the momentum equations. We find that 



dcf 



h ^' solution to 
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;iS°.6^ = sp-ffvuJoG 



cos(T) 

— sin(T) 

/5cos(z/T) 

(3 



sTiojqG 



ox^ 



ox^ 



P\x^ 



(3.20) 



/ 



We note that ox'^H^ph^ = = oi"|: (s^^g^'^). 

We write x" = ox" + /i li" and solve |: (Mx") = eF°'^x^ -hP''f'£ {^^^^3b'^)■ Substituting 
for P"^ and rearranging terms yields -^ (Mi" + UTj'^ph^] = eF°'pX^ — hx^'x^-^ f S^^6^) = 
eF^^i;^ - /i oi" ox''£ {^^^pb^) = eF^^i^ since qx'^^ (s«^6^) =0. For a = i = 1, 2 we see 
that 



(m — s/icjo) i* = eF\x^ . 



(3.21) 



This is independent of g, but dependent on the magnitude of the spin, and leads to a shift 
in the cyclotron frequency. This shift is in the other direction for a particle with spin down. 



eBz 



m—jhujo 



— jl^. This effect may 



For a particle with s = | one sees that 'yuJcydotron 

be of interest for the case of plasma dynamics in critical magnetic fields of the order of 

10^^ — lO^^Gauss, which may exist near pulsars. 



Continuing with the analysis we find that i*^ = 7 = constant and 



X 



j2 hsuo 
m 



— /3^ z/cos(z/T) + constant 



dr \ m 



(3.22) 



where for simplicity we have set the integration constant equal to zero. The charge will 
radiate in virtue of this oscillation along the axis of the applied magnetic field. The oscillation 
exists because the four-spin is orthogonal to the four-velocity, and the spin precesses. 

IV. CONCLUSION 



The theory predicts a contribution to the electric dipole radiation of the particle with 



frequency co = uJo{^ ~ 1) and power on the order of |cA;^ |p|^ 



6c3 



/3^e^) due to the 



' 2mcl 
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0{h) oscillations of the orbit along the direction of the applied static magnetic field. (There 
are of course other 0{h) sources of electric dipole radiation that we are not considering 
here. One such source is proportional to fida^j^-) This suggests a new way to measure 
the gyromagnetic ratio of the electron or proton. One can measure | — 1 by measuring the 
frequency of the electric dipole radiation due to the 0{h) oscillations of the orbit along the 
direction of the applied static magnetic field. 

Radiation reaction terms have not been included in the analysis. However, radia- 
tion damping tends to polarize electrons in high energy accelerators. These effects are 
observable and must be included in any complete description of relativistic dynamics of 
light particles. The radiation damping force is usually assumed to arise from the self- 
field of the accelerated charged particle via the Lorentz force. Dirac has shown that 
the finite (on the world line of the particle) contribution to the self-field is given by 
F^f^ = \ (F",f - F^"j^) = f (v'^v^ - t;^i)") = f |: (v'^i)'^ - w^w"), where t;" = i". Adding 
this field to the externally applied electromagnetic field F"^ in Equation [ |2.12| ] yields a gen- 



eralization of the Abraham-Lorentz equation. Such equations are well known to be plagued 
with "self accelerated" runaway solutions. One simple way around this is to replace w" with 
-pF'^f^, since an external F^"^ is ultimately responsible for ?)" 7^ 0. This yields a second 
order ODE for x'^ (coupled to the first order ODE for ■?/;). 
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